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Исследованы задачи распределения неоднородных работ между неоднородными исполнителями с учетом 
затрат на выполнение работ и переключение между ними. Предложены генетические алгоритмы распреде- 
ления работ, которые могут использоваться при решении различных практических задач как в автоматизи- 
рованных процессах, так и в социальных системах. 
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Введение. Задачи распределения работ между несколькими исполнителями возникают в различ- 
ных областях науки, техники и социальных сферах, причем характер работ и правила их распре- 
деления могут быть различными в зависимости от рассматриваемой прикладной области. 

Как правило, задачи распределения работ относятся к сложным многоэкстремальным за- 

дачам, эффективные методы решения которых известны только для ряда частных случаев. Обыч- 
но для получения точных решений применяют метод ветвей и границ, метод динамического про- 
граммирования или эвристические методы [1, 2]. Эффективность метода ветвей и границ в суще- 
ственной степени зависит от «точности» граничных оценок подмножеств решений. Известные эв- 
ристические методы и метод динамического программирования применимы только для опреде- 
ленного класса или даже вида задач и неэффективны или неприменимы при решении наиболее 
общей и наиболее приближенной к практике задачи — распределения работ в случае, когда за- 
траты зависят от последовательности выполнения работ. Поэтому актуальной является разработ- 
ка эффективных методов решения задач распределения неоднородных работ между несколькими 
неоднородными исполнителями. 
Постановка задачи. При распределении работ между исполнителями необходимо минимизиро- 
вать затраты на их выполнение. В понятие затраты включаются потери времени, расход электро- 
энергии, износ оборудования и т.д. В целях перехода к однокритериальной задаче оптимизации 
предлагается свести все учитываемые затраты к одному параметру. (Способ перехода от много- 
критериальной задачи к однокритериальной в данной статье не рассматривается). 


Пусть имеется п (: 1) неоднородных исполнителей и т ( 7 =1т) неоднородных работ, 


представленных в виде множеств: А = {х....,",} и И ={м,...*„} соответственно. 


2 т 

Особенностью рассматриваемой задачи является то, что процесс переключения исполни- 
теля с работы на работу сопровождается затратами, которые зависят как от исполнителя, так и от 
последовательности работ. 

Необходимо получить последовательность выполнения заданных работ так, чтобы при за- 
данных свойствах исполнителей и работ минимизировать затраты на выполнение работ с учетом 
затрат на переключение. 

Задания поступают на вход системы динамически, т.е. в любой момент времени { на вход 
системы может поступить задание. Таким образом, необходимо перераспределять работы в ре- 
альном масштабе времени, которые на момент времени 2 не выполнены или не выполняются. То 
есть работы находящиеся на выполнении не подлежат замещению. 

Исполнитель характеризуется затратами на выполнение данного вида работы и затратами 
на переключение с одного вида работ на другой. Исполнитель может быть не способен выполнить 
какие-либо работы. 
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Затраты на выполнение исполнителями определенных работ представляются в виде мат- 
рицы 7 размерности пхт, в которой элемент 2, показывает затраты на выполнение /-й работы 


1-м исполнителем. А затраты на переключение исполнителя с работы на работу - в виде матриц 
Р® размерности тхт (и). 


Таким образом, в результате решения задачи необходимо получить такое распределение 
работ между исполнителями, чтобы обеспечить выполнение следующего критерия: 


п т Ш т т 
- (ОО. 
Е=>, а, > руиу > ша, (1) 
2 = 1 у 
где с, - элемент матрицы С размерности пхт , принимающий значение 1 или 0 и показываю- 


щий, выделена ли исполнителю 7 для выполнения работа /; иФ — элемент матрицы И“ 
размерности тхи ‚ принимающий значение 1 или 0 и показывающий, планируется ли для 


1-го исполнителя переключение с работы / на работу /. 
Первое слагаемое выражения (1) определяет затраты на выполнение всех планируемых 


работ, а второе - затраты на переключение исполнителей с работы на работу. 

Задачу предлагается решать с использованием генетического алгоритма (ГА), который 
представляет собой простую модель эволюции в природе, реализованную в виде компьютерной 
программы [3]. В нем используются как аналог механизма генетического наследования, так и 
аналог естественного отбора. При этом сохраняется биологическая терминология в упрощенном 
виде. 

Алгоритм решения задачи распределения работ зависит от количества исполнителей 

(один или несколько) и от того, учитываются ли затраты на переключение с работы на работу. 
В связи с этим предлагается провести декомпозицию задачи и получить алгоритм формирования 
последовательности выполнения работ для одного исполнителя и алгоритм распределения работ 
между исполнителями без учета затрат на переключение, а затем, путем объединения указанных 
алгоритмов, записать общий алгоритм. 
Алгоритм формирования последовательности выполнения работ для одного исполни- 
теля (Алгоритм 1). В случае одного исполнителя задача сводится к минимизации стоимости вы- 
полнения всех работ и аналогична задаче коммивояжера без необходимости возвращения в ис- 
ходный пункт. При этом необходимо получить такую последовательность выполнения работ, что- 
бы выполнялся критерий: 


т т 
Е=У Ури, > п. 


ИГУ 
Целевая функция (функция фитнеса) определяется как ф=-Ё и наилучшей является 


хромосома, имеющая наибольшее значение функции фитнеса. 

Как известно, решение с помощью ГА сводится к последовательному выполнению сле- 
дующих шагов: 

— выбор символьной модели объекта в виде закодированных хромосом; 

— генерация начальной популяции; 

— модификация хромосом в целях улучшения популяции; 

— выбор хромосомы из популяции с максимальным значением функции фитнеса. 

Способ кодирования определяется спецификой задачи. При выборе символьной модели 
необходимо обеспечить возможность кодирования в хромосоме любой точки из пространства по- 
иска. Невыполнение этого условия может привести к невозможности найти решение поставлен- 
ной задачи. 
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Хромосома отражает последовательность выполнения работ исполнителем. Как указано 
выше, исполнитель должен выполнить т работ. Каждая хромосома С’ представляет собой конка- 


тенацию ряда подкомпонентов, называемых генами С = (2,2.... 2) ‚ гГДе т - длина хромосомы. 


Число генов в хромосоме равно числу работ, поступивших на выполнение. Ген с номером 1 опре- 
деляет номер работы, которую необходимо выполнить 1-й в последовательности. 

Поскольку одна работа должна быть выполнена только один раз, то каждая работа может 
входить в хромосому единожды. 

Алгоритм получения начальной популяции заключается в следующем. Начальная популя- 
ция включает м хромосом. Первым геном 1-й хромосомы является 1-я работа. Каждым последую- 


щим /-м (7=2,т) геном является номер работы с минимальными затратами на переключение 


между текущей и последующей в соответствии с жадным алгоритмом. 

Из полученных хромосом выбирается хромосома с максимальным значением функции 
фитнеса, т.е. хромосома, которая кодирует последовательность работ с минимальными затратами. 

Затем осуществляется модификация выбранной хромосомы с помощью операции мутации. 

В данном алгоритме не используется операция кроссовера, а только операция мутации 
над единственной выбранной хромосомой. 

Можно предложить несколько операций мутации. 

1. Перестановка двух произвольно выбранных генов. На рисунке приведен пример мута- 
ции путем перестановки генов 2 и 4. 


1 2 3 4 5 6 
ЕЕ 
а) 

1 2 3 4 5 6 
вены 

5) 


Пример операции мутации, заключающейся в перестановке двух произвольных генов: 
а- хромосома до мутации; б- хромосома после мутации 


2. Циклический сдвиг на одну позицию генов, находящихся между двумя случайно вы- 
бранными. 

3. Перестановка последовательности генов в другое место хромосомы. 

4. Замена генов на обратную последовательность. 

5. Перестановка последовательности генов в другое место хромосомы в обратной после- 
довательности. 

При модификации выбранной хромосомы предлагается повторять все пять мутаций с рав- 
ной вероятностью до тех пор, пока не будет получено решение с удовлетворительными затратами 
или в течение некоторого количества мутаций не будет найдено решение с меньшими затратами. 

Для увеличения вероятности правильного выбора переставляемых генов и уменьшения 
времени работы алгоритма введем понятие тесносвязанных работ (ТСР). Это такие работы г и }, 
для которых затраты на переключение с работы 1 на работу 7 меньше или равны затратам пере- 
ключения с работы г на любую другую, кроме той, что предшествует работе г, или с любой другой 
на работу у. Соответственно, если эти условия не выполняются, то работы называются нетеснос- 
вязанными (НТСР). 
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Таким образом, разрыв хромосомы предлагается проводить таким образом, чтобы «не 
разрывались» тесносвязанные работы. Перестановка может осуществляться в любую позицию. 
После перестановки необходимо заново определить точки разрыва. 

Если полученная после модификации последовательность работ, описанная хромосомой, 
имеет меньшие затраты, то она принимается вместо исходной, если нет - то отбрасывается. 

В качестве примера проведем распределение работ, затраты на переключение между ко- 
торыми представлены матрицей: 

о 10 020 30 40 50 150 140 


10 0 1020 30 4 150 5 
20 0 0 0 2003090 30 
30 20 0 0 002 98 9 
40 30 20 10 0 100 90 10 
50 40 30 20 0 0 10 80 
150 150 60 50 90 10 0 10 
100 90 80 100 110 120 10 0 


Необходимо выполнить каждую работу единожды. Поэтому затраты на выполнение каж- 
дой работы неважны. 

Воспользовавшись алгоритмом получения начальной популяции, определим следующий 
набор хромосом: 





12345687; 21345687; 32145687; 43215687; 
54321687; 65432187; 78321456; 87432156. 
В качестве начального приближения выберем первую хромосому 12345687, поскольку она 
кодирует последовательность работ с минимальными затратами на переключение, равными 140. 
Обозначим НТСР знаком «-». В результате получим хромосому 123456-87 с двумя группа- 
ми ТСР. 
В результате выполнения операции мутации получена хромосома 123-87—456 путем при- 
менения модификатора 3. Затраты на переключение для этой хромосомы равны 130. 
Предложенный метод решения имеет существенное преимущество перед ГА без введения 
понятия ТСР и полным перебором. Ясно, что количество операций будет зависеть как от количе- 
ства генов (количества распределяемых работ), так и от количества ТСР. 
Сравнительная характеристика полученного алгоритма с классическим ГА приведена в 
таблице. 
Сравнительная характеристика алгоритмов 



































Количество Количество Время Время Относительная погрешность по- 
генов НТСР работы ГА, мс работы ГА с ТСР, мс лученного решения ГА с ТСР, % 
8 3 63 12 0 
8 6 63 25 0 
24 6 270 31 0 
24 12 270 56 0 
24 16 270 89 1,1 
32 6 634 34 0 
32 12 634 64 1,4 
32 18 634 93 2,4 
100 18 3102 127 2,2 
100 42 3102 1136 2,8 























Из последнего столбца таблицы видно, что предложенный генетический алгоритм с ТСР 
может находить решение, отличное от точного с достаточно малой погрешностью. 
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Алгоритм распределения работ между исполнителями без учета затрат на переключе- 
ние (Алгоритм 2). Для распределения работ между несколькими исполнителями также можно 
использовать ГА. В этом случае необходимо получить такое распределение работ, чтобы выпол- 
нялся критерий: 
п т 
Е, = а: —>шш. 
1=Г = 
Целевая функция (функция фитнеса) определяется как ф=-Ё, и наилучшей является 


хромосома, имеющая наибольшее значение функции фитнеса. 

Предлагается следующий алгоритм генерации начальной популяции при количестве ис- 
полнителей равном и. Начальная популяция будет включать хромосомы в количестве, кратном 
количеству исполнителей. Кодировка хромосомы имеет вид: в г-м гене находится номер исполни- 
теля, который выполняет 1-ю работу. Гены /-Й хромосомы содержат номер исполнителя 
7 то4 т, где шо4 - операция вычисления остатка от целочисленного деления / на тя. В случае, 


если исполнитель / не может выполнить работу 1, то в гене г во всех хромосомах популяции необ- 
ходимо заменить число / на номер другого исполнителя. Таким образом, можно учесть предпоч- 
тение передачи работы какому-то исполнителю. 

Поскольку один исполнитель может выполнить несколько работ, то номер одного и того 
же исполнителя может присутствовать в разных генах. 

При этом важным является правильный выбор размера популяций. Размер популяции мо- 
жет повлиять как на скорость сходимости, так и на качество найденного решения. Правиль- 
ный выбор размера популяций может быть осуществлен только экспериментальным путем. 
В проведенных экспериментах размер популяции для 8 исполнителей и 32 работ принимался в 
количестве 80. 

Проиллюстрируем алгоритм получения начальной популяции на примере. Пусть имеется 3 
исполнителя и 6 работ. 

Примем размер начальной популяции равным шести хромосомам. В случае, если любой 
исполнитель может выполнить любую работу, то начальная популяция имеет вид: 

111111; 222222; 333333; 111111; 222222; 333333. 

Для случая, если исполнитель 3 не может выполнить работу 2, передадим работу 2 испол- 
нителям 1 и 2 с равной вероятностью. В связи с этим в хромосомах 3 и 6 работа 2 передана ис- 
полнителям 1 и 2 соответственно, что обеспечивает равную вероятность выполнения работы 2 
этими исполнителями. Таким образом, начальная популяция имеет вид: 

111111; 222222; 313333; 111111; 222222; 323333. 

Генерация начальной популяции является подготовительным этапом к процессу субопти- 
мизации, который заключается в последовательном применении к полученной популяции генети- 
ческих операторов согласно ГА. 

В данном алгоритме используется модель ГА Уитли (0. МПЖЕеу). При скрещивании произ- 
вольно выбранных двух хромосом применяется одно- или двухточечный кроссинговер с равной 
вероятностью. Местоположение точек разрыва выбирается случайным образом. В результате 
скрещивания получается один потомок, к которому применяется операция мутации. 

Мутация заключается в замене в каком-либо гене номера исполнителя на другой, выбран- 
ный случайно. Количество мутаций предлагается принимать не более 5% от количества работ. 
Так, при 32 работах были выполнены мутации в 1-2 генах. 

При выполнении Г скрещиваний полученная хромосома занимает в популяции во время 
первых [/2 скрещиваний место родительской хромосомы с минимальным значением функции 
фитнеса, а на последних [//2 скрещиваниях - место хромосомы с минимальным значением функ- 
ции фитнеса во всей популяции. 
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Алгоритм заканчивает свою работу при достижении удовлетворительного решения или 
при выполнении указанного количества скрещиваний. При проведении вычислительных экспери- 
ментов при 32 работах, 8 исполнителях и размере популяции равном 32 количество скрещиваний 
принималось равным 750, что позволяло получать решение, близкое к оптимальному. 

Алгоритм распределения работ между исполнителями с учетом затрат на переключе- 
ние (Алгоритм 3). Этот ГА представляет собой комбинацию двух вышеприведенных алгоритмов 
и позволяет минимизировать критерий (1). Алгоритм имеет вид. 

1. Генерация начальной популяции в соответствии с Алгоритмом 2. 

2. Определение последовательности работ для каждой хромосомы популяции в соответ- 
ствии с Алгоритмом 1 и вычисление значения функции фитнеса. 

3. Сохранение максимального значения функции фитнеса ф 


тах" 

4. Отбор двух родителей и получение потомка путем скрещивания и мутации в соответст- 
вии с Алгоритмом 2. 

5. Вычисление значения функции фитнеса ф для полученного потомка. 


6. Если количество полученных потомков менее [/2, то замена потомком родителя с ми- 


нимальным значением функции фитнеса, иначе — замена потомком наименее приспособленной 
хромосомы в популяции. 


7. Если ф>ф 
8. Генерация начальной популяции для каждого исполнителя в соответствии с Алгорит- 


то перейти на шаг 8, иначе -— перейти на шаг 4. 


глах / 


мом 1. 
9. Получение последовательности работ в соответствии с Алгоритмом 1. 
10. Вычисление значения функции фитнеса ‘у с учетом последовательности работ. 


11. Если \у>ф„„, то принять фи, =. 


тах 7 

12. Если достигнуто удовлетворительное решение или выполнено указанное количество 
скрещиваний, то перейти на шаг 13, иначе — перейти на шаг 4. 

13. Конец алгоритма. 

После получения хромосомы, описывающей распределение работ между исполнителями 
без учета затрат на переключения, для каждого исполнителя выполняется алгоритм формирова- 
ния последовательности работ для одного исполнителя с учетом затрат на переключение. Как 
следует из алгоритма, в целях экономии времени алгоритм формирования последовательности 
работ для каждого исполнителя (шаги 8-12) выполняется в том случае, если полученная хромо- 
сома имеет значение функции фитнеса больше сохраненного максимального значения ф 


тах" 
Только в этом случае есть возможность получить распределение работ между исполнителями с 
затратами меньше ранее найденного решения. 

Для уменьшения времени работы алгоритма также используется понятие ТСР. Однако это 
понятие в данном случае принимает следующий вид. Тесносвязанными являются работы ги /, для 
которых затраты на переключение с работы 1 на работу / меньше или равны затратам переключе- 
ния с работы { на любую другую этого исполнителя, кроме той, что предшествует работе г, или с 
любой другой работы этого исполнителя на работу /. 

Проведенные эксперименты показали, что относительная погрешность найденного реше- 
ния по сравнению с оптимальным для восьми исполнителей и 32 работ при количестве скрещива- 
ний равном 320 не превышает 7%. При этом время работы алгоритма составило около 9 с на со- 
временном ПК. 

Заключение. Предложенные алгоритмы позволяют осуществить распределение работ между 
исполнителями и сформировать последовательность выполнения работ с учетом затрат на пере- 
ключение. Проведенные вычислительные эксперименты подтверждают возможность использова- 
ния полученных алгоритмов при решении различных практических задач, связанных с распреде- 


лением работ, как в автоматизированных процессах, так и в социальных системах. Последующие 
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доработки связаны с проведением исследований и получением рекомендаций по выбору парамет- 
ров генетического алгоритма, обеспечивающих уменьшение вычислительных затрат, сходимость 
алгоритма и нахождение приемлемого субоптимального решения. 
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